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SUMMARY 


The  problem  of  calculating  airplane  stability  parameters  from  the 
aircraft  response  to  an  arbitrary  disturbance  is  considered.  To  calcu¬ 
late  the  coefficients  of  the  linear  differential  equation  which  describes 
the  airplane  transient  response,  application  is  made  of  a  classical 
least- squares  curve  fitting  method. 

It  is  shown  that  the  method  is  applicable,  although  somewhat, 
cumbersome,  when  the  input  is  an  arbitrary  function  of  time.  Certain 
inputs  are  demonstrated  to  lead  to  simplification  in  the  application  of 
the  method.  Examples  are  given  illustrating  the  means  of  using  the 
method  and  showing  its  practicability. 

Finally,  an  appendix,  in  which  Prony’ s  method  and  a  generalization 
thereof  appears,  is  presented. 


INTRODUCTION 


The  determination  of  the  stability  and  control  parameters  of  a 
dynamical  system  from  its  measured  response  has  been  a  subject  of 
increasing  importance  and  interest,  both  from  the  standpoint  of  basic 
aerodynamic  research  and  automatic  stabilization  of  airplanes.  This 
problem  is  essentially  one  of  curve  fitting  and  may  be  stated  as  follows: 
Given  the  time  history  of  an  airplane  response  to  a  known  transient  dis¬ 
turbance,  and  assuming  a  certain  form  of  linear  differential  equation 
with  constant  coefficients  describing  the  relation  between  response  and 
disturbance,  required  to  find  the  coefficients  of  this  differential 
equation  such  that  the  sum  of  the  squares  of  the  differences  between  the 
given  response  and  the  one  corresponding  to  the  differential  equation  is 
a  minimum. 
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A  number  of  methods  for  the  solution  of  this  problem  have  been 
advanced  by  others  and  are  noted  in  reference  1.  However,  these  methods 
do  not  apply  the  least  squares  principle  in  the  correct  sense  described 
above . 

The  method  presented  herein  involves  obtaining  the  solution  of  the 
differential  equation  in  a  form  amenable  to  rigorous  application  of  the 
method  of  least  squares  and  such  that  a  completely  arbitrary  input  may 
be  treated.  A  number  of  means  are  available  for  the  fitting  of  a 
function  in  this  form  by  least  squares.  The  solution  by  means  of  a 
Taylor’s  series  expansion  seemingly  gives  good  results  and  is  used 
throughout  this  paper. 

Considerable  simplification  of  the  method  is  possible  when  free 
oscillation  data  are  available  or  when  the  input  can  be  exactly  or  at 
least  closely  represented  by  a  function  the  analytical  form  of  which  is 
known.  Solutions  for  a  pulse,  a  step,  or  a  ramp  input  are  given  in 
detail  or  are  indicated.  The  general  solution  for  an  arbitrary  input  is 
also  included. 

Selected  examples  of  different  inputs  applied  to  the  same  physical 
system  (airplane)  are  presented,  which  illustrate  the  different  means  of 
applying  the  method. 

The  report  has  been  so  organized  that  the  engineer  who  is  not 
interested  in  the  derivation  of  the  formulas  used  may  read  the  first 
section  of  the  report  entitled  "Statement  of  the  Problem”  and  proceed 
from  there  directly  to  the  section  on  examples. 


METHOD  OF  ANALYSIS 


Statement  of  the  Problem 


r  >  Consider  a  quantity  qm(t)  which  has  been  measured  at  a  set 
i  tip  i  =  0,  1,  .  .  .  ,  V~1  of  V  values  of  t.  Suppose  the  initial 
conditions  on  qm,  obtained  from  the  experimental  data,  are 


where 


n  is  determined  below. 
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A  physical  interpretation  of  these  quantities  may  he  had  hy  con¬ 
sidering  the  data  qm(t)  as  a  time  history  of  the  pitching  velocity  of 
an  airplane  in  response  to  an  elevator  deflection  F(t). 

Consider  next  the  differential  equation 

Po(D)q(  t)  =  Pi(D)F(  t)  (2) 


where  D  is  the  operator  d/dt,  PQ  =  Dn  +  an_1  D0-1  +  .  .  .  +  aQ  and 
Pi  =  CjhD01  +  C^i  I)111”1  +  .  .  .  +  CD  are  two  polynomials  in  D, 

and  F(t)  is  a  known  forcing  function.  Consider  the  set  of  all  possible 
solutions  of  this  differential  equation  obtained  by  varying  the  constants 
ai  and  C  j .  Let  qc(t)  be  that  solution  of  equation  (2)  for  which 


M  = 


v-l 

[q(ti)  -  qm(ti)]2 


i=0 


(3) 


a 


is  a  minimum,  and  subject  to  the  same  initial  conditions  as  qm 


qc(°)  =  *o  ] 


J 


It  is  then  desired  to  find  those  values  of  the  a^  and  the  Cj  which 
correspond  to  qc.  These  values  will  be  unique,  since  to  each  set  of 
these  constants  there  corresponds  one  and  only  one  solution  of  the  differ¬ 
ential  equation  if  the  initial  conditions  are  determined. 

In  most  practical  problems,  m  is  less  than  n,  and  this  assump¬ 
tion  will  be  made  throughout  this  paper.  The  extension  of  the  method 
described  herein  to  the  case  where  m  is  equal  to  or  greater  than  n 
should  be  clear. 
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Let  Xi,  \2,  .  .  Xn  "be  the  zeros  of  the  polynomial  Pg(d).  It 
is  now  assumed  that  A-p  f  A.j,  i  ^  j*  The  extension  of*  "the  following 

solution  to  the  case  where  not  all  the  \j[  are  distinct  can  be  found 
in  any  textbook  on  elementary  differential  equations  and  needs  not  be 
included  here  (reference  2).  It  is  well  known  that  if 


Pof(D)  =  ^  Po(D) 


the  solution  of  equation  (2)  is 


q(t) 


n 


.  +  m 

0  VUj) 


i 


J-l 


e  J  F(  t  )  dT 


w 


where  the  constants  Aj  are  functions  of  the  initial  conditions  q(o). 


and  of  the  constants  Ci,  C2>  .  .  .  ,  Cm.  The  method  described  in  this 
paper  consists  of  fitting  q(t),  by  a  least  squares  procedure,  to  a 
function  of  the  form  (4).  Several  methods  for  fitting  such  a  function 
may  be  found  in  the  literature.  In  this  report  the  classical  device  of 
linearization  and  iteration  by  means  of  a  Taylor fs  series  (reference  3> 
p.  214)  is  used.  Therefore,  a  detailed  study  of  this  linearization 
follows,  but  it  must  be  understood  that  it  is  in  no  way  essential  to  the 
method.  Fitting  by  steepest  descents  (reference  4),  for  example,  could 
be  used  in  place  of  the  Taylor* s  series  iteration.  The  mechanics  of  this 
iteration  are  cumbersome  in  the  general  case.  They  will,  therefore,  be 
described  for  particular  cases,  thus,  it  is  hoped,  making  the  method 
clear  for  any  special  case  which  may  arise. 


Formulation  of  the  Method  When  Free 
Oscillation  Data  Are  Available 


The  method.—  Suppose  there  exists  a  T  such  that  F(t)=0  for  all 
t  >  T.  Then,  for  t  >  T,  the  expression 


Ai 


pi(M 
o'Uj) 


F(t  )dr 


occurring  in  equation  (4)  is  constant.  Call  this  constant  B^. 

J 
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Then  for  t  >  T 


n 

q(t)  =  ^  Bj  eXjt  (5) 

j=l 


4 


By  some  means  (see  e.g.,  reference  l),  first  approximations  to  the 
constants  occurring  in  equation  (6)  are  found.  Let  B2^°^ , 

.  .  .  ,  Bn(°^  Xi(°),  .  .  .  ,  “be  the  approximations  to 

Bi,  BS,  .  .  .  ,  Bn,  Xi,  \2,  .  .  .  ,  \n>  respectively.  The  function 
q(t)  is  now  expanded  by  means  of  a  Taylor  series  about  the  "point” 

\bx(°\  .  .  .  ,  Bn(°),  X^o),  .  .  .  ,  Xn(o) 


with  all  terms  of  order  higher  than  the  first  omitted.  From  equation  (5), 


_c)q 

dBi 


X-t  t  dq 
e  1  ,  — 
c)Xi 


Bi  t 


i  =  1, 
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Denote  the  values  of  q,  at  the  "point"  .  .  .  , 

“  (.rVis-V01.  ^ia.V0) 


•  •  •  >  X] 


(o) 


By  q’ 


'n  1 

Then  by  Taylor’s  theorem,  omitting  all  higher  order  terms. 


,  respectively. 


where 


-t(i 


(o) 


n 


i=l 


SiY  -  *  I 

i=l 

Aq  =  q  -  q(°) 

ABi  =  Bi  -  Bj/0) 
AXi  =  Xi  -  Xi<°) 


{tf  - 


(6) 


i=l,  .  .  .  ,  n 
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Therefore,  from  equation  (6) 


M  = 


=  \  e  1  ABi  +  Bj_'  °'  t  AXi 

i=l 


The  desired  equations  which  lead  to  the  minimization  of 
V— 1 

12 


(7) 


9c(*i)  “  Mti) 


i=0 


may  now  he  found  as  follows:  Let 
^m  =  9m  -  9(o) 


Then 


I 


^fl(ti)  -  A3m(ti) 


i=0 


The  equations  for  the  increments  then  become 


8m 

d(ABi) 


dM 

^AX*) 


>(8) 


i  1,  •  •  .  y  n  J 

Since  the  expression  (7)  for  Aq  is  linear  in  (ABi)  and  (AXi), 
the  equations  (8)  will  be  precisely  those  found  when  fitting  redundant* 
data  to  a  linear  function  by  least  squares  (reference  3,  pp.  209-211). 
Equations  (8)  are  solved  for  the  increments  ABi  and  AXi,  and  the 
resulting  values  are  added  to  Bi^0)  and  X^0),  respectively,  to 
give  the  second  approximation. 


The  iteration  process  by  means  of  Taylor’s  expansions  is  now 
repeated  until  two  successive  iterations  give  the  same  values  for  the 
parameters  to  the  desired  number  of  significant  figures. 
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It  is  of  course  true  that  in  many  cases,  especially  those  in  which 
very  good  data  are  available,  no  such  Taylor  series  expansion  is  neces¬ 
sary,  Prony’s  method  (appendix  A)  giving  sufficiently  accurate  values  of 
the  parameters.  Examples  have  been  found,  however,  where  Prony!  s  method 
did  not  lead  to  sufficiently  good  values,  and  the  Taylor  series  was 
essential. 

Final  calculation  of  the  coefficients  of  P0(D)  and  Pi(D).—  All 
that  now  remains  for  the  problem  to  be  completely  solved  is  the  computa¬ 
tion  of  the  coefficients  of  P0(l))  and  Pi(D)  from  the  calculated 
values  of  Bi  *  and  \±. 

Since  P0(D)  =  Dn  +  a^^  D3^1  +  .  .  .  +  a0,  and  since  Xi,  .  .  .,  Xn 
are  the  zeros  of  this  polynomial. 


n 


(reference  5>  P-  29). 

The  constants  Bj_  were  defined  as  follows: 


Bi  =  Ai  + 


PiUi) 

Po’Ui) 


e  ^-*-T  F(t)cIt 


where  T  is  that  value  of  t  where  free  oscillation  begins.  It  is 
well  known  that  the  A^  may  be  solved  for  as  functions  of  the  constants 
Ci,  .  .  .  ,  Cm.  The  integrals  occurring  in  the  definition  of  Bj_  can 
now  be  found  graphically  since  X±  and  F(t)  are  known.  There  are, 
therefore,  n  equations  in  the  m-1  unknowns  CQ,  •  •  .  ,  Cm  (the 
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Tiie  Method  for  a  Second— Order  System 


A  case  of  great  practical  interest  in  aerodynamics  is  that  of  the 

dSn  71'^  to„the  Pr!sent  s"‘1”  ^  -ill  ie  Seated  S 

detail .  A  second-order  system  is  one  in  which  equation  (l)  becomes 


(D2  -  a2D  -  a0)  q(t)  =  (CXD  -  CQ)  F(t) 


An  example  of  such  a  second-order  system  which  occurs  in  aero' 
dynamics  my  be  obtained  by  letting  ,  be  the  response^  p?tc5S 
velocity  to  an  elevator  input  F(t).  In  this  case  (reference  1), 


=  iL 

mV0  Iy  Iy 


=  -i*L 

mVrt 


MS  Lfr  Le. 

C0  =  -X- ± - “■  x 


Iy  mV0  I, 


mV, 


CX  =  ^  -  IS.  x  ^ 

Iy  mVo  Iy 


where 


m  mass  of  airplane 

V0  trim  velocity  of  airplane  * 

Iy  pitching  moment  of  inertia 
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L  lift  force 

M  pitching  moment 

5  elevator  input  F(tt) 

a  angle  of  attack 


a 


da, 

dt 


ht 


MS 


dl 

dS 

di 

da 

dM 


Ma 

Mdt 

Mq 


dM 

da 

dM 

dd 

dM 

3q 


The  general  solution  of  (D 2  +  axD  +  a0)q  = 
obtained  from  equation  (4)  which  becomes,  if  Xi 
of  x2  +  axx  +  aQ  =  0, 


q  =  AxeXlt  +  A2eXat  +  eXlt 

ki“A,2 


C3A2+C0 


e-X2T  F(  t  )  dT 


Differentiating,  an  expression  for 
letting  t=0,  it  may  be  seen  that 


Ai  +  A2  -  <l(  0) 


(CxD  +  C0)F  may  be 
and  are  the  roots 


e  ^lT  F(t  )dr  + 


may  be  founds  and 


AiXx  +  A2\2  —  4(0)  —  C  iF(  0 ) 
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or,  using  equation  (la). 


Therefore, 


qo-qoXar-CJ'o 
Xx- "X2 

q0-io^-CJ0 

XgfXx 


40-<io^2-Ci^'0+(c1^i+Co)  Io  e  ^lT  5’(t)^t  exxt  + 

Xi““ A.  2 


$  — X2T 


qQ-^ qWq^^  Iq  e  ex2t 


Pulse  input.—  When  a  pulse— type  input  has  been  applied  to  a  second- 
order  system,  that  is,  when  there  is  a  T  such  that  F(t)=0  for  all 
t  >  T, 


Bx  -  Ai  + 


B2  =  A2  + 


Ci^i+0o 

Xx“A/2 


X2  ^1 


e  ^iT  f(t)cLt 


e-X2T  F(T)dT 


\  (10) 


using  the  same  notation  as  before.  Then,  for  t  >  T, 


=  Bx  eXlt  +  B2 


Suppose  now,  as  is  generally  the  case,  that  Xx  and  X2  are  con¬ 
jugate  complex  numbers.  The  form  previously  used  for  the  description  of 
the  general  method  may  still  be  followed  in  this  case;  however,  simplifi¬ 
cations  occur  if  this  notation  is  abandoned.  If  Xi  and  X2  are  con¬ 
jugate  complex,  Bx  stud  b2  must  be  so  also,  since  q  is  real.  Let 
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Then 


\x  =  l  +  V i 
Bx  =  p  +  p*i 

\2  =  Z  -  1*1,  B2  =  p  -  p'i 


and 

q  =  2elt  (p  cos  l' t  -  p*  sin  Z't)  (12) 


It  is  assumed  that  at  this  point,  first  approximations  ZQ,  Z0* , 
PQ,  P0'  to  the  constants  Z,  V ,  p,  p',  respectively,  have  heen 
obtained.  (See  appendix  A  for  one  method  for  finding  such  approxima¬ 
tions  . ) 

The  formulas  needed  for  the  Taylor's  series  iteration  will  now  he 
derived.  From  equation  (12), 


V— 1 

■I 

i=0 

Taking  the  derivatives  of  M  with  respect  to  AZ,  AZ' ,  Zip,  Zip'  and 
setting  them  equal  to  zero  will  result  in  the  equations,  the  solution  of 
which  will  give  those  values  of  the  increments  which  minimize  M. 

These  equations  thus  become 
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As  before.  It  is  nov  assumed  that  this  iteration  is  repeated  until 
two  successive  iterations  give  the  same  values  for  the  parameters  to  the 
desired  number  of  significant  figures „ 
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The  parameters  Z,  V ,  0  and  0*  have  thus  "been  determined.  The 
solution  to  our  problem  will  be  complete  if  we  can,  from  these  values, 
calculate  aQ,  ax,  Co,  Cj..  From  equations  ( 9) , 


ax  —  +’ ^2  “  2Z 


^2  =  ^1^2  =  22  +  Zt2 


A  means  of  finding  Ci  and  Co  must  nov  be  found.  In  order  to 
accomplish  this,  consider 

T  T 

J  e~^iT  F(r)dr  =  J  e“*T  cos  Vt  F(t)cLt 


Let 


Let 


Then 


T 

— i  j  e~ sin  Z 1  t  F(t  )dT 
Jo 

cr  =  /  e~^T  cos  Z't  F(r)dT 

Jo 

r T 

a'  =  /  e  sin  Z’t  F(r)dT 

e~^lT  F(t  )dr  =  a  —  a'  ± 


and  since  X2  is  the  conjugate  of 


e  ^2T  F(t  )dr 


a  +  cr*  ^ 


Recalling  that  Bi  =  p  +  P'i,  it  follows  from  equation  ( 4a)  that 
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-  (  )Ci-crtC0+Z,q0 

21' 

b»  =  Zqo~^o~[lcr+z,(T,-F(0)  3ci-qCn 

2V 


These  equations  may  "be  solved  for  CD  and  Ci. 

The  method  when  the  forcing  function  continues  throughout  the 
motion.-  It  is  now  supposed  that  in  the  given  data,  there  is  no T  such 
that  the  forcing  function  F(t)=0  for  all  t  >  T.  In  this  case,  the 
coefficient  of  e  J  in  the  expression  for  q(t)  given  by  equation  (4) 
is  no  longer  constant.  Prony’ s  method  cannot,  therefore,  be  used  for 
the  first  approximation.  The  simplification  which  occurs  in  the  Taylor' s 
expansion  in  the  previous  case  due  to  the  fact  that  q(t)  is  a  sum  of 
exponentials  with  constant  coefficients  is  also  not  present  here. 

I 

There  are  cases,  however,  even  when  free  oscillation  does  not  occur 
where  the  present  method  can  be  applied  with  little  difficulty.  The 
best  examples  of  such  cases  are  those  in  which  F(t)  is  known  to  have  a 
certain  analytical  form.  As  an  example,  the  two  degrees  of  freedom 
system  (D2  +  bD  +  k)  q(t)  =  (CjD  +  CQ)  F(t)  may  be  considered,  where 

F(t)  =  t  >  0#  so_ca-*--*-e<^  "step"  function. 


For  this  input,  equation  (4a)  becomes 


q  _  ^-i4q— ^lkagp+Cp  exxt  +  ^aio— ^lkagp+Cq  gx2t  +  Co 

^2)  k2(X2—A.i)  A.lA.2 


(4a) 


where,  since  dq/dt  is  discontinuous  at  t=0  if  F(t)  is  a  step, 

3*0  =  TT  >  >  0+  indicating  that  t  is  to  approach  zero  through 

positive  values  only.  The  apparent  anomaly  that  q  is  independent  of 
Ci  is  resolved  in  appendix  B,  where  it  is  shown  that 


4(0+)  =  q(0— )  +  CiF(0+)  =  4(0-)  +  Cx 
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First  approximations  to  the  constants 


Xi(\i— ^2) 


Xxq(  0—)— Xi\2{1o+CiXi+Cq 
Xx(  X1—X2) 


B  _  ^4(  0+)-"XiX2qo+Co 

X2(  X2— Xi) 


Xgq(  0— )~X1X2qQ+C1X2+CQ 
Xa(  X2—X1) 


which  occur  in  equation  (4b)  can  be  found  by  a  simple  extension  to 
Prony1  s  method.  (See  appendix  A.)  A  Taylor fs  expansion  may  then  be 
applied  to  obtain  a  closer  approximation. 

A  very  important  input,  closely  related  to  the  step  is  the  so-called 
"pseudo— step.11  In  any  physical  case,  an  exact  step  can  never  be  obtained 
for  an  input.  The  input  will  always  have  a  certain  finite  slope  near 
t=0;  there  may  also  be  a  certain  amount  of  overshoot  or,  reciprocally, 
the  input  may  undershoot  its  steady  state  value.  Any  input  which  rapidly 
(but  not  instantaneously)  attains  a  constant  nonzero  value  will  be  called 
a  pseudo-step.  Suppose  then  that  F(t)  is  such  a  function.  Then  there 
is  a  value  T  such  that  F(t)  is  a  constant  for  all  t  >  T.  Suppose 
this  constant  is  c.  Let  this  function  be  applied  to  the  second— order 
system  (D2  +  ax  +  a0)  q(t)  =  (CXD  +  CQ)  F(t).  The  response  in  q  is 
then  given  by  equation  (4a).  It  is  here  assumed  for  simplicity  that  the 
initial  conditions  q(0)  and  q(0)  are  both  zero;  F(o)  is  also  zero. 
Then 


q(t) 


CjAi+cq  e\xt  rt 

J  n 


e~XlT  F(  t  ) d t  +  — 

X2~-Xi 


e  f(t  )dT 


For  t  >  T, 


q(t) 


CiX,i+C0  Xit 

.  e 

Xx — X2 


r  T 

F  e  ^lT  F(t)(1t  + 

F  e  ^1T  F(T)d  r 

Jo 

T 

+ 


CxXg+Cp  x2t 

X2”XX 


e  ^2  t  f(  t  )  d  t  + 


e  X2T  F(  T )  d  T 
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ip  — \  .  q- 

fQ  e  F(t)<1t  is  constant  and  since  F(t)  is  constant  for  t  >  T,  nnfl 
jt  e  Xit  ^  t)^t  may  -be  evaluated. 


c  CiXi+Cq  _  p  c  Ci^+Co  ^ 

Xi  xr-x2  “  Ql>  TPkT  “  2' 


CiX'i+Cq 

Xl“—\2 


^  _W  ,  V  CiXa+Co  pT 

e  xiT  F(r)dr  =  Hx,  — — - —  /  e~x2T  F(r)dT  =  H2 

0  XjHtl  J0 


Then  Gi,  G2,  Hi,  H2  are  constants,  and 


q(t)  =  Hi  eXit  +  Gi  +  H2  e^  +  G2  £eM  'fc-T)-l 


=(Hl*J^)  ellt+(Ia  +  x?) 


e^2  —  (Gi  +  G2) 


Thus  q(t)  is  the  sum  of  two  exponentials  plus  a  constant.  First 
approximations  may  now  he  found  hy  the  extension  to  Prony' s  method  given 
in  appendix  A,  and  Taylor's  series  expansions  may  he  used  to  improve 
these  values. 


Another  simple  example  which  may  he  considered  is  that  where 

F( t)  =  {  ^  ®  .  In  this  case  the  response  of  the  second— order  system 

would  he* x>  t  “  0 


q  =  7ieXlt  +  72e%st  -  {7iXi+7s^z)  t  -  ( 7i+7z ) 


q(0)  =  q(0)  =  0 


vhere 


_  Ci\i+Co  _  CiXgtCp 

^i2(  Xi  X>2 )  ^22(  X2r“Xi) 
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A  first  approximation  may  be  found  by  several  means ,  the  derivative 
method  (see  example  3)  being  only  one  of  these.  Again,  Taylor’s  series 
may  be  applied  to  improve  these  values. 

In  the  most  general  case,  where  F(t)  is  only  known  graphically  or 
tabularly,  the  method  can  still  be  applied;  there  is  a  certain  amount 
of  added  labor  involved,  but  in  case  such  an  example  should  occur,  the 
method  for  this  general  case  is  outlined  below.  It  is  to  be  noted  that 
this  general  case  is  better  suited  to  computation  with  a  high-speed 
electronic  computer  than  it  is  to  manual  computation.  However,  it  may 
certainly  be  applied  manually. 

The  entire  method  proceeds  from  equation  (4).  It  is  assumed  that 
the  given  differential  equation  has  been  solved  for  q  in  this  form. 

Step  1.—  By  some  means,  a  first  approximation  to  the  parameters 
must  be  f bund  (several  such  methods  may  be  found  in  reference  l)  . 

Step  2.—  The  function  q(t)  is  expanded  in  a  Taylor’s  series,  all 
teams  of  order  higher  than  the  first  being  omitted.  Herein  lies  most  of 

.  t  j  «j- 

the  computation,  since  J  Q  e  F(t )dT  must  be  found  graphically  as  a 
function  of  t. 

As  an  example,  consider  again  the  second— order  case  (D2  +  axD  +  ao) 
q(t)  =  (CiD  +  C0)  F(t),  where  F(t)  is  given  graphically.  It  is 
assumed  for  simplicity  that 

f0  =  q(o)  =  q(o)  =  0 

Then,  using  the  same  notation  as  before,  Ax  =  A2  =  0,  and  q(t)  is 
simply  the  particular  integral: 


q(t)  =  ellt  f t  e^lTF(T)4T  +  Sjg£S>  e*®*  ^  F<T>dT 

Jq 


The  computations  (especially  when  it  comes  to  taking  derivatives  for  use 
in  the  Taylor  series)  are  greatly  simplified  by  the  substitution 
_  C  C2.A,2+Cn 

=  - x  B2  =  - -  .  The  problem  is  then  restated  as  follows: 

To  find  the  best  fit  for  q(t)  to  a  function  of  the  form 


q(t)  =  Bi  eXlt  f  e  XlT  F(t )dT  +  B2  eXst  e~X2T  F(t )dT 


o 
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1 ✓ 

Minimization  of  M  =  (qc— qm)  with  respect  to  3,  3f,  Z,  Z1,  (where 

3  +  3*i  =  Bi,  3  —  3 1  i  Z  +  Z’i  =  Xx,  Z  —  Z'i  =  Xz)  is  equivalent 

to  minimization  with  respect  to  Ci,  C0,  a.±,  and  ao,  which  is  the 
desired  minimization.  The  method  now  proceeds  as  before,  with  the  first 
approximations  being  found  and  with  the  application  of  a  Taylor1  s  series . 

Another  means  which  might  be  employed  in  this  general  case  would  be 
to  first  fit  F(t)  to  some  suitable  function,  and  then  find  q(t)  as 
an  explicit  function  of  t  as  was  done  above  for  the  step  and  the  ramp. 


EXAMPLES 


The  engineer  who  is  interested  primarily  in  applications  may  read 
these  examples  immediately  after  reading  the  first  section  of  the  report 
entitled  "Statement  of  the  Problem." 


Example  I  ~  Pulse  Input 


Consider  the  input  shown  in  figure  1(a)  applied  to  the  secondr-order 
system 


(D2  +  arf)  +  ao)q(t)  =  (CJ)  +  Co)F(t)  (16) 


The  response  is  given  in  figure  1(b) .  The  methods  and  formulas  described 
previously  can  be  used  directly  for  this  example.  If  t  >  0.4, 


q  =  Bx 


Xxt 


+  B; 


\zt 


(5a) 


or,  if  Bj,  =  0  +  0'i,  Xi 


+  2*  i,  and  t  >  0.4, 


q 

2 


e^  (0  cos  2 •  t  —  0’  sin  2't) 


(12) 
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It  is  now  assumed  that  by  some  means  first  approximation  Z0,  Z0' , 
Po,  Pc',  to  Z,  V  y  0’ ,  respectively,  have  been  found.1  Prony' s  method 
(see  appendix  A)  gives  Z0  =  —  O.9I,  Z0'  =  7*02,  PD  =  O.368I, 

P0*  =  -  2.6775. 

From  equation  (12) 


8Z 


t  eU  (p  cos  Z't  -  p»  sin  Z't) 


t  eZ^  (P  sin  Z't  +  P’  cos  Z't) 


—  f—  =  eZ^  cos  Z't 

dp  \2/ 


UD 


elt  sin  Z't 


Prony' s  method,  for  example,  may  be  used  here .  Prony' s  method  is 
actually  applied  to  this  example  in  appendix  A.  An  ordinary  period 
and  damping  analysis  (described  below),  which  is  shorter  than  Prony' s 
method,  may  also  be  used.  The  afore-mentioned  period  and  damping  anal¬ 
ysis  proceeds  as  follows:  The  minima  of  q  lie  on  a  curve  q  =  K  e^ 
K  =  constant.  Let  qi  and  q2  be  two  successive  minima  of  q.  Then 

(-—S\=  e^T.  t  =  period  of  a  (which  mav  be  found  from  ficmre  1  _ 

w 


herefore. 


t  =  period  of 
I  =  A  lnfS*N 


(which  may  be  found  from  figure  1(b)). 
As  for  Z ' ,  Z '  is  equal  to  —  .  Then 


q  =  2  e (p  cos  Z't  —  p’  sin  Z't) 


Fitting  at  any  two  points  tx  and  t2 

(cos  Z'tx)  p  -  (sin  Z'ti)  p'  =  -  e-Ztl  q(tx) 

2 

(cos  Z't2)  P  —  (sin  Z't2)  P'  =  —  e— Z^2  q(t2) 

which  may  be  solved  for  p  and  p' .  This  analysis  gives  Z  £-0.9, 

Z'  ~7»0.  Letting  ti  =  1.0,  t2  =  2.0,  it  may  be  found  that  p  ~0.:+, 
p’  «2 .6. 
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Therefore,  by  Taylor's  theorem,  omitting  all  higher  order  terms, 

t(p  cos  Z't  -  p'  sin  Z't)  AZ  -  t  (p  sin  Z't  +  p»  cos  Z't)  AZ’  + 

(cos  Z't)  AP  -  (sin  Z't)  AP'  =  —  e~lt  -  (p  cos  Z't  -  p*  sin  Z't) 

Referring  to  table  I  vhere  circled  numbers  refer  to  columns, 

©  AZ  -  @  AZ'  +  ©Ap  -  ©  Ap'  =  © 

Minimizing  M  leads  to  the  equations  (see  equations  (l4) ) 

AZ  £©2-AZ'  Z  @  @  +  AP  E  ®  @-AP»  E  @  @  =  £  ©  © 

-AZ  £  ©  @+AZ'  E  ©2-  AP  £©©  +  AP'  2  ©  ©  =-  2  @  © 
Al  2@  ©-AZ'E  @  @+  APE  ©2-Ap'  2  @©  =  Z@  © 

-AZ  Z©  ©+  AZ'  Z  ©  @-AP£©©  +  AP'£  ©2=-  2  ©  © 
These  equations  give 

AZ  =  -  0.0100  AZ'  =  0.0052 

Ap  =  -  0.0118  AP'  =  -  0.0320 
or 

Z  =  -  0.9200'  Z'  =  7.0252 
P  =  0.3563  P'  =  -  2.7095 


A  second  iteration  was  now  applied,  and  it  was  found  that  the  increments 
were  zero  to  four  decimals. 
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Finally,  for  computing  a1(  aQ,  Ci  and  C0 


ax  =  -  2Z  =  1.84 

a0  =  Z2  +  Z ' 2  =  50.1998 

The  functions  e— ^  F(t)  cos  Z't  and  e—i^F(t)  sin  Z’t  are  now  tabu¬ 
lated  and  plotted  (table  II  and  fig.  2).  These  two  functions  must  now 
be  integrated  to  find  a  and  cr*  .  It  is  to  be  noted  that  in  this 
example  F(t)  is  exactly  equal  to  t  for  0  <  t  <  0.2,  and  to 
(0.4  —  t)  for  0.2  <  t  <  0.4,  and  so  the  integration  may  be  done 
analytically.  The  integration  may  be  done  graphically  by  means  of  a 
planimeter,  however,  when  F(t)  is  not  so  clearly  a  perfect  triangular 
pulse,  as  is  usually  the  case.  A  planimeter  was  actually  used  in  this 
example  to  give  a  =  0.00493  and  o'  =  0.0405.  Equations  (15)  are  now 
set  up: 

0.280  Cx  +  0.00493  C0  =  38.0690 
0.07189  Ci  -  o.o4o5  c0  =  5.0061 

which  give  Ci  ='  134.0  and  CQ  =  114.2,  to  four  significant  figures. 
The  true  values  of  the  parameters  in  the  preceding  example  were 


ax  =  1.84 
aQ  =  50.2 
Cx  =  134.0 
C0  =  114.4 


Example  II  —  Step  Input 


Consider  a  "step"  input 
described  by  equation  (l6). 


F(t)  =  ^  t  >  0  aPPlied  to  the  system 
The  response  is  given  by  equation  (4b). 


From  the  response  shown  in  figure  3  It  may  be  seen  that  q(0)  =  0, 
q(0— )  =  0.  Then  equation  (4b)  becomes 
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q(t) 


C3.X3.+C0  x,  t  Ci\2+Co  \„t  cQ 

— - - -  e  x  +  — - - -  e  +  - 

X2)  ^1^2 


* 


Let 


3  -  CjAi+Cq 

B  =  ^Xg+Cp 


Then,  using  the  same  notation  as  before,  with  \x  =  2  +  Zfi,  A.2  =  l  —  VI, 
Xi  +  X2  ~  “  &l,  ^1X2  ”  ao> 


Co  Ciao+CoZ 
2a0  22* a0 


Cp  +  Ciap+CgZ  ^ 
2a0  22’ aQ 


Let  Bx  =  p 
Therefore, 


+  Pfi. 

B3  =  - 


Then  *>  ■  -  a?  - 

2g  .  Then  |  =  ezt  ($ 


Cia0+C02 
22' a0  * 


But  B3 


Co 

ao 


cos  2't  —  p'  sin  2* t)  —  p  . 


First  approximations  may  now  Be  found  by  the  extension  to  Prony' s 
method  described  in  appendix  A  or  by  a  generalization  of  the  ordinary 
period  and  damping  analysis.  For  simplicity,  this  period  and  damping 
method  was  used  to  give  2Q  =  -  0.91,  20'  =  7.0,  pQ  =  -  1.11,  p0*  =  -9.7  . 


tr  |  oj 
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which  in  turn  give 
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l  =  —  0.9200 
Z'  =  7.0242 
P  =  -  1.1392 
P*  =  -  9.3890 

A  second  Taylor's  expansion  was  then  applied  (see  table  IV ) ,  giving 

l  =  -  0.9201 
l •  =  7.0258 
p  =  —  1.1400 
P*  =  -  9.3875 


Another  iteration  would  lead  to  increments  which  are  zero  to  four 
decimals . 

For  the  final  calculation  of  the  coefficients. 


ax  =  -  (A.1+X2)  =  -  2 Z  =  1.8402 
aD  =  X1X2  =  J2*'2  =  50.2084 
CQ  =  2kp  =  114.4715 


c x  =  -  =  134.0072 

k 


It  is  to  he  noted  that  a  second  iteration  such  as  was  applied  above 
was  hardly  important,  since  the  results  of  the  first  Taylor's  series 
give 
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ax  =  1.84 
ao  =  50.19 
C0  =  114.34 

Cx  =  134.00 

the  true  values  1)61115,  as  before, 

ax  =  1.84 
aQ  =  50.2 
CQ  =  114.4 
Ci  =  134.0 

Example  III  —  Input  Requiring  the  General  Method 

Consider  finally  the  input  shown  in  figure  4(a)  applied  to  the 
system 

(D2  +  aiD  +  aD)q(t)  =  (CJ)  +  C0)F(t)  (1 6) 


Such  an  input  might  occur,  for  example,  in  a  stabilized  airplane  where 
the  pitching  velocity  q  is  fed  back  to  the  elevator  to  change  the 
input. 

First  approximations  by  the  derivative  method.-  The  so-called 
"derivative  method"  will  here  be  applied  in  order to  find  first  approxi¬ 
mations  to  the  desired  constants.  From  figure  4,  F  and  q  are  found 
graphically  as  functions  of  t.  Then  q  is  plotted  (fig.  5),  and  from 
this,  q  is  found  and  tabulated.  (See  table  V.)  Rewriting  equation  (16), 


(q)ai  +  (q)ao  -  (F)Ci  -  (F)C0  =  -  (q) 
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and  referring  to  table  V, 

(f)ai  +  (3)ao  -  (5)Ci  —  0CO  =  —  (j) 


A  least-squares  analysis  is  now  applied  to  this  equation,  giving 


ax  =  1.84 
a0  =  50.19 

Ci  =  133.89 
C0  =  114.91 


The  high  accuracy  of  these  first  approximations  is  due,  of  course,  to  the 
excellence  of  the  data  in  the  example.  The  leading  objection  to  the 
derivative  method  obviously  is  the  necessity  of  finding  the  derivatives 
graphically,  which  in  many  cases  leads  to  errors  so  gross  as  to  make  the 
values  of  the  parameters  found  in  this  manner  entirely  valueless. 


The  Taylor's  series  iteration.-  Since  F(t)  is  continuous  at  t=0, 
the  constants  q(0),  4(C),  F(0)  of  equation  (4a)  may  be  found  as  usual, 
by  inspection  of  figures  4  and  5.  Since  q(0)  =  q(0)  =  F(0)  =  0,  from 
equation  (4a), 


q 


C  1A.1+C0 
Xl“A/2 


p(T)dT  +  ^ 

^2-^-1  J0 


ex2T 


F(T)dT 


or,  letting  Bi 


CiXi+Cq 

Xi—X2 


B2 


CiX2+Cq 
X.2-X-X  ' 


q  =  Bi 


F(r)dr  +  B2 


F(r)dT 


Using  the  same  notation  as  before,  with  \i  =  l  +  Z'i,  Bi  =  0  +  p'i, 
/J*  e“ f(t)(1t  =  cr  —  a' i,  it  follows  that  X2  =  2  —  Z’i, 

B2  =  p  -  P'i,  /q  e— ^2T  F(r)dT  =0  +  0*1,  and  that. 
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V 


% 


q  =  2  eJt[(0o  +  0'd')  cos  Z't  +  (per*  —  p * cr)  sin  Z't]  (17) 


It  will  also  "be  noted  that  0  =  ^  ,  0'  =  -  .  Letting  zero  sub- 

2  2Z' 

scripts  denote  the  values  of  the  indicated  parameters  at  the  first 
approximat ions , 


ZQ  =  “  0.92 


Va  =  7.02 


P0  =  66.945 


P’o  =  -  0.59 

From  equation  (17) 

~  =■  2  eZt  -f  [  (ptx+p*  cr*  )t  -  (pp+p*p»)]  cos  Z't  + 
oZ  L 

[(0o'-0'o)t  +  (P'p-Pp')]  sin  Z't  j- 
-^3.  =  2  elt\  [  (per'— p* cr) t  +  (p'p-0p*)]  cos  Z't  - 

bv  l 

[  (Pcr+P'cr'  )t  -  (0p+0'p')]  sin  Z'tj- 


vhere 


—  =  2  (cr  cos  ZTt  +  cr1  sin  Z’t) 

dp 

~  =  2  eU  (o'  cos  Z't  -  o  sin  Z’t) 

dp' 


r 


p=  P  t  r(t)  e  ^  cos  Z't  dt 
o 

p'  =  P  t  F(t)  e— ^  sin  Z't  dt 
'-'o 


(18) 
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The  function  F(t)  is  now  multiplied  by  the  four  quantities 
e”~^0^  cos  lf0t,  e— sin  V  0t,  t  e  cos  Z*0t,  t  e  sin  Z’0t, 


point  by  point,  and  the  resulting  products  are  plotted  against  t 
(fig.  6).  These  four  curves  are  integrated  as  functions  of  t,  giving, 
respectively,  cr,  crf ,  p,  and  p1  (table  Vi). 


The  quantities 


qp  (dq/dl)0  (bgJbV  )Q 

^0^  * 


(Bq/Sp)0 

,  exot  ' 


are  then  computed  from  equation  (18). 


and 


Then,  using  a  Taylor1  s  series  and  table  VI, 


©A  l  +  @AV  +  ©A3  +  ©A3’  =  © 
Minimizing  M  with  respect  to  A l,  Al' ,  A3*  A3'  gives 


A 1  =  -  0.00 

AV  =  0.01 


A3  =  0.17 
A3*  =  -  0.05 

Due  to  the  inherent  inaccuracy  of  a  planimeter,  only  two  decimal  places 
were  preserved  here. 

Thus 


l  =  -  0.92 

v  =  7.03 

3  =  67.03 
3»  =  -  0.64 
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The  desired  parameters  were  finally  computed  to  he: 

a^_  =  1.84 

a0  =  50.28 

Ci  =  134.06 
C0  =  114.69 

CONCLUDING  REMARKS 

A  method  has  been  described  by  which  the  coefficients  of  the  differ¬ 
ential  equation 


dHq  dn  1q  dq 

3  *  ■>«  3?=* +  •  •  •  * ai  i +  ■»» 


d1^  dm~1F 

m  dtm  m"1  it131-1 


^  dF  _  _ 
+  Ci  —  +  CqF 
dt 


can  be  calculated  from  the  graphical  knowledge  of  q(t)  and  F(t).  It 
is  noticed  that  the  method  may  become  somewhat  cumbersome  if  the  input 
F(t)  is  not  of  certain  types.  The  input  which  allows  the  method  to  be 
applied  most  easily  is  one  which  goes  to  zero  quickly,  giving  free  oscil¬ 
lation  data.  A  simplification  also  occurs  when  the  input  F(t)  is 
known  accurately  to  be  a  function  of  a  certain  type,  such  as  a  step,  a 
ramp  F(t)  =  t  or,  perhaps,  something  of  the  form  £  Aie^ 

The  method  has  as  its  primary  advantages  first,  the  fact  that  the 
correct  quantity  (the  sum  of  the  squares  of  the  differences  between  the 
calculated  and  the  measured  quantities)  is  minimized.  The  ordinary 
simple  equations  of  least  squares  may  be  used  in  this  minimization,  since 
the  equations  of  condition  (the  redundant  or  inconsistent  equations  to 
be  solved  by  least  squares)  satisfy  the  restrictions  under  which  the 
least  squares  solution  is  normally  derived.  That  is,  the  equations  are 
linear,  and  only  the  right-hand  sides  are  subject  to  error.  The  so- 
called  "derivative  method"  used  in  example  III  minimizes 
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..  •  <=: 

Z(q  +  1d(5  +  kq  ~  CjF  —  CqF)  ,  a  quantity  with  little,  if  any  meaning. 
Other  methods,  Prony's,  the  Laplace  and  Fourier  methods,  are  subject  to 
the  same  objection. 

Another  advantage  of  the  method  is  that  the  data  are  analyzed 
directly  in  the  time  plane,  thus  eliminating  the  possibility  of  the 
introduction  of  errors  due  to  graphical  integration  (in  the  case  of 
Fourier  and  Laplace  transform  methods)  or  differentiation  (derivative 
method) .  Where  the  inputs  give  free  oscillation  data  there  is  a  certain 
saving  of  time  due  to  this  directness,  no  preliminary  steps  having  to  be 
taken  before  the  method  can  be  applied. 

For  the  usual  type  of  problem  which  is  found,  those  with  the  pulse- 
or  step-type  inputs,  the  method  is  at  least  as  rapid  as  any  other  which 
gives  comparable  results. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Jan.  ±7,  1951. 
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APPENDIX  A 

„  PRONY'S  METHOD 


A  problem  which  occurs  quite  often  is  that  of  fitting  a  function 
q(t)  to  a  sum  of  exponentials.  One  method  used  to  fit  a  function  of 
this  form  is  known  as  Prony's  method  and  is  described  below.  (See  also 
references  1  and  3.) 


It  is  first  assumed  that  the  data  are  given  in  tabular  form  with 
equal  intervals  of  the  argument  t.  Suppose  measurements  of  q(t)  are 
taken  at  instants  t0,  t*,  .  .  .  ,  ty— v  >  2n,  where  n  is  the  number 
of  exponentials  to  which  q  is  to  be  fitted.  Then  it  is  assumed  that 
tk  -  t0  +  k(At),  At  =  constant,  k  =  0,  1,  .  .  .  ,  v-1.  By  a  proper  shift 
of  the  time  axis,  t0  may  be  taken  equal  to  zero.  Let  qk  =  q(tk)  = 

q(kAt).  Fundamental  to  Prony's  method  is  the  following  theorem:  If 

n 


k  =  0,  1, 
tion 


v-l. 


i=l 


Bi  e^k 


(19) 


then  q  satisfies  the  linear  difference  equa— 


<lk+n  +  Qnr-i  ak+IML  +  %-z  <3k+n_,  +  .  .  .  +  Qoqfc  =  0  (20) 


where  Qq,  .  .  .  ,  Qn_1  are  constants  such  that  the  roots  of  the  equa¬ 
tion 


*n  +  Q^i  xD_1  +  •  •  •  Qq  =  0 

Xi(At) 

are  e  >i-l,  #..,n.  The  proof  is  as  follows: 


i-L 

9k+n  +  Qn-i  ^k+n— 1  "*■•••+  =  ^  Bi  eM.(k+n)At  + 

i=l 

“  n 

Qn-i  £  Bi  e^**-1^*  +  .  .  .  +  Qo  V  Bi  e^ 


(21) 


i=l 


B^  e 


XiAt 


’WAt)°  *  ^  ♦ . . .  ♦ 


i=l 

n-i 


Qo 


i=l 
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Certainly  there  exist  n  constants  Q_.,  ....  Q-.  . ,  such  that 
XiAt  X2At  XnAt  ^  ^i_1 

e  >  e  y  •  •  •  )  e  are  the  roots  of  equation  (21) .  Choosing 

these  constants  in  this  way  then  makes  the  bracketed  expression  above 
and,  therefore,  qjj-+n  +  On-!  <lk+n-i  +  .  .  .  +  Q0  vanish.  This  completes 
the  proof . 

Prony* s  method  consists  of  writing  down  the  equations 

^n  +  ^n— 1  ^n— 1  +  •  •  •  +  Qo^o  =  ^ 

^n+i  +  1  ^n  +  •  •  •  +  ^0^1  =  ® 

^n+2  +  ^n— 1  ^n+i  +  •  •  •  +  ^0^2  =  0 


and  solving  them  by  least  squares2  for  Qi,  i  =  0,  .  .  .  ,  n-1.  From 
the  normal  equations  obtained  from  the  leasb-squares  process, 

Q0>  •  •  .  9  Qxb-i  can  found.  Then  the  roots  of  the  equation 

x11  +  Qn-x  x11""1  +  .  .  .  +  Q0  =  0 

are  calculated,  giving  e^i(At)  ^<3^  therefore, 

2Herein  lies  one  of  the  objections  to  Prony1  s  method,  since 
v-n 

^  ^k+n  +  Qn-i  ^k+n~i  +  •  •  •  +  Q03k) 

L _ I 

k=0 

rather  than 

'?  (■.-£»  -'‘"X 

k=0  i=l  1 

is  minimized,  the  correct  minimization  procedure  leading,  in  this 
case,  to  a  forbidding  amount  of  calculation. 
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The  coefficient  e*'1*'  of  Bi  in  equation  (19)  can  he  tabulated 
for  each,  i,  since  is  known.  Therefore,  can  he  found  hy  a 

second  least— squares  procedure.  This  completes  Prony1 s  method. 

The  objection  described  in  the  preceding  footnote  can  be  overcome 
by  considering  the  and  found  by  Prony1  s  method  not  as  the  best 

possible  values  of  these  parameters,  but  only  as  a  first  approximation 
thereof.  A  better  approximation  can  then  be  found  by  means  of  a  Taylor1  s 
series.  This  method  is  described  in  the  body  of  this  report. 


Extension  to  Prony1 s  Method 


Suppose  now  that  it  is  required  to  find  some  way  to  fit  a  function 
to  a  sum  of  exponentials  plus  an  (unknown)  constant.  The  solution  will 
be  presented  herein  for  the  case  in  which  the  exponentials  are  two  in 
number,  but  the  generalization  to  a  greater  number  of  unknowns  will  be 
evident . 

To  be  specific,  a  quantity  q  must  be  fitted  to  a  function  of  the 

form 


q  =  Bi  e^  +  B2  eX2t  +  B3  (22) 

q  is  first  fitted  to  a  function  of  the  form 

q  =  Ba  eXlt  +  B2  eX2t  +  B3  eXst  (23) 

and  the  condition  that  \3  =  0  is  put  in  later.  Prony' s  method  is 
applied  first  to  equation  (2) .  As  before  q  satisfies  the  difference 
equation 

^k+3  +  Q2  qk+2  +  ^k+i  +  Qo^k  =  0  ( 2 ^  » 

where  the  constants  Q*,  i=0,  1,  2,  3,  are  such  that  the  equation 

x3  +  Q2  x2  +  Qx  x  +  Qo  =  0  (25) 

has  as  roots.  But  X3  =  0.  Therefore,  equation  (4)  has 

unity  as  a  root,  and  1  +  Q2  +  Qx  +  Q0  =  0.  Eliminating  Q0  (say) 
between  this  and  equation  (24) 
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2  “  ^•'k)  ^2  +  (^k+i  ^k)  ^1  (^k  ^k+3^  {%&) 

is  obtained.  Thio  equation  is  now  solved  by  least  squares  for  Qx 
and  Q2.  The  exponentials  e^-^t)  and  are  now  found  to  be 

the  roots  of  x2  +  (Q2  +  l)  x  +  (Qx  +  Q2  +  1)  =0.  Pronyfs  method  nov 
proceeds  as  before. 

An  evident  extension  of  this  method  may  be  used  to  simplify  the 
problem  of  fitting  a  sum  of  exponentials  when  one  or  more  of  the  expo¬ 
nents  are  already  knovn. 


Example  of  Prony*  s  Method 

Suppose  the  data  found  in  column  (^3)  ,  table  VII  is  to  be  fitted  to 
tvo  exponentials  (see  example  I  of  the  body  of  this  report)  •  In  this 
case,  equations  (20 )  become 

@+0Q1+(3)Qo=O 

Solving  by  least  squares , 

^Q1Z@2+Qo2@  ©  =  -  2  ©  (D 

.Qi2®  ©+Qo2:(3)2=-£®  © 


which  give 


Q0  =  0.8320 
Qx  =  -  1.3922 


Equation  (21)  thus  "becomes 


x2  -  1.3922  x  +  0.8320  =  0 


or 
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=  eXl^At^  =  O.696I  +  O.5895  i 


Therefore, 


Xi  =  —  In  (O.696I  +  O.5895  1) 

At 


=  10 


2n  „/(0.696l)2+(0.5895)2  +  1  arc  tan  ^ *5^95 

O.6961 


=  -  0.91  +  7.02  1 


X2  =  -  0.91  -  7.02  i 


Since  Xx  and  X2  are  complex  conjugate,  so  are  Bx  and  B2.  Let 
=  P  +  P'i.  Then  q  =  2  elt  (p  cos  2't  —  p*  sin  2't).  Referring  to 
table  VII,  this  becomes  (j)  P  —  (8)  P'  =  (12)  .  Solving  as  before  by- 
least  squares, 

P  =  O.368I 


3‘  =  -  2.6775 
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APPENDIX  B 

DISCONTINUOUS  INPUTS 


Consider  the  differential  equation 


a-a 

dtn 


+ 


dn~1q 

&a~1  dt13-1  + 


+  ai  «  +  S°q  =  °m  ^ 


+ 


+  C 


dF 

dt 


+  C  oF 


(2) 


A  fundamental  question  which  often  arises  is  the  following.  Suppose 
F(t)  or  one  of  its  derivatives  is  discontinuous  at  a  point.  At  this 
point  all  higher  derivatives  fail  to  exist.  What  is  the  meaning,  if  any, 
of  the  differential  equation  (2)  at  this  point? 


An  instance  of  this  problem  occurs  in  example  II  of  the  present 
report.  There,  F(t)  is  a  step:  F(t)  ={  J'  £>  0  *  At  *  =  0 
dF/dt  and  all  higher  derivatives  of  F  do  not  exist. 

Throughout  the  following  discussion,  let  x(0+)  =  lim  x(t), 

t  — >  0+ 

where  x(t)  is  any  function,  and  t  -5>  0+  indicate  that  t  is  to 
approach  zero  through  positive  values  only.  Consider  again  the  differ 
ential  equation 

q  +  &i  q  +  a0  q  =  CXF  +  C0F  (l6) 


where  dots  denote  differentiation  with  respect  to  t.  Integrate  equa¬ 
tion  (1 6)  between  the  limits  -e  and  t  (e>0)  to  obtain 


[q(t)  -  q (— e )  ]  +  ax[q(t)  -  q(-e)] 


q(t)dt 


=  Ci  [  F(t)— F(— e)  ] 


+  c< 


F(t)dt 


r  4 

It  is  now  assumed  that  J_eq(t)dt  is  continuous,  (in  all  examples  used 
in  this  report,  q(t)  represents  the  pitching  velocity  of  the  airplane, 
making  ^  q ( t ) dt  the  angle  of  pitch  which,  from  physical  considera¬ 

tions,  may  be  seen  to  be  continuous.) 
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F(t)dt  = 


F(t)dt 


since 


t  <  0 
t  >  0 


Letting  €  — >  0, 

[4 ( t ) — 4(0—)]  +  ax  [q(t)~ q(O-)  ]  +  k[0 (t)-O(O- )  ]  «  CiF(t)  +  C0t 


vhere 


0(t) 


Integrating  once  more  with  respect  to  t,  this  time  from  —  G  to  +G, 


[<i(g)— q(— e)— 2eq(0-)  ]  +  ax  [0(e)— d(— e)— 2eq(0-)  ]  + 


0(t)dt  -  2e@(0— ) 


=  Ci  G 


Again  letting  e — >0,  it  is  seen  that  q(0+)  —  q(O-)  =  0,  since 
0(t)  was  assumed  continuous.  Thais  q(t)  is  continuous  at  zero.  Going 
hack  to  equation  (19)  and  letting  t — >0,  q(0+)  —  q(O-)  =  CiF(0+)  =  Ci. 
Thus,  q(t)  is  discontinuous  at  zero,  and  the  difference  between  the 
right-hand  and  left-hand  limits  of  q(t)  as  t  approaches  zero  is  Ci. 
Thus,  writing  equation  (4a)  in  terms  of  q(0— )  rather  than  q(0+), 

,  \  kiq(0— )— XxX2q(0)+Ciki+Co  x^t  X2(l(0— )~ ^•i^2(l(0)+CiA.2+Co  X2t  , 

9(t) - - - e  * - UtPZT) - e  + 


and  q(t)  may  "be  seen  after  all  to  le  dependent  on  C*. 
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TABLE  I.-  THE  TAYLOR’S  SERIES  ITERATION  APPLIED  TO  EXAMPLE 
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TABLE  II.-  FURTHER  CALCULATIONS  NECESSARY  FOR  EXAMPLE  I 


2 

3 

4 

5 

6 

7 

8 

9 

10 

t 

0.92t 

e— .92t 

7.0252t 

COS  (5) 

sin© 

F(t) 

©@® 

G3S!! 

0 

0 

1.0000 

0 

1. 0000 

0 

1 

0 

0 

0 

.05 

.046 

.9550 

.3513 

.9389 

.3441 

.05 

.04483 

.0164 

.1 

.092 

.9121 

.7025 

.7632 

.6461 

.1 

.06961 

.0589 

.15 

.138 

.8711 

1.0538 

.4943 

.8693 

.15 

.06469 

.1136 

.2 

.184 

.8319 

1.4050 

.1650 

.9863 

.2 

.02745 

.1641 

•25 

.230 

.7945 

1.7563 

-.1844 

.9828 

.15 

-.02198 

.1171 

.3 

.276 

.7588 

2.1075 

-.5113 

.8594 

.1 

-.03880 

.0652 

.35 

.322 

.7247 

2.4588 

-.7758 

.6310 

.05 

-.02811 

.0229 

.4 

.368 

.6921 

2.8100 

-.9455 

.3256 

0 

0 

0 
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TABLE  III.-  TAYLOR'S  SERIES  ITERATION  APPLIED  TO  EXAMPLE  II 
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TABLE  V.-  TABULATION  OF  F(t),  F(t),  q(t),  <£( t ) j,  AND  q(t)  FOR  EXAMPLE  III 


■9 1 

M3M 

mmmm 

1 

t 

F 

F 

q. 

4 

q. 

I 

0 

0 

7.03 

0 

0 

941 

2 

.1 

.5893 

4.35 

3.940 

68.4 

327 

3 

.2 

.8206 

.21 

10.939 

57-1 

-532 

k 

.3 

.6521 

-3.33 

12.974 

-23.9 

-978 

5 

.4 

.2251 

-4.81 

5.847 

-113.9 

-702 

6 

.5 

-.2290 

-3.92 

-7.837 

-146.1 

no 

7 

.6 

-.5062 

-1.46 

-20.407 

-91.3 

938 

8 

•  7 

-. 5l4l 

I.23 

-24.019 

24.6 

1266 

9 

.8 

-.2947 

2.92 

-15.576 

137.9 

886 

10 

•  9 

.0175 

3.05 

1.300 

184.8 

5 

11 

1.0 

.2694 

I.81 

18.214 

138.3 

-895 

12 

1.1 

.3606 

-.01 

26.577 

21.5 

-1334 

13 

1.2 

.2779 

-1.53 

22.181 

-105.4 

-1092 

l4 

1.3 

.0869 

-2.11 

7.319 

-178.7 

-312 

15 

1.4 

-.  1103 

-1.67 

-10.555 

-163.6 

594 

16 

1.5 

-.2258 

-.57 

-22.775 

-71.0 

1171 

17 

1.6 

-.2226 

.60 

-23.785 

50.6 

1155 

18 

1-7 

1221 

1.31 

-13.711 

l4l.4 

589 

19 

1.8 

.0152 

1.32 

2.053 

160.2 

-219 

20 

1.9 

.1228 

.75 

15.764 

103.1 

-866 

21 

2.0 

.1582 

-.05 

21.207 

2.5 

-1058 

22 

2.1 

.1182 

-.70 

16.504 

-91.2 

-741 

23 

2.2 

.0329 

-.93 

4. 663 

-134.9 

-107 

24 

2.3 

-.0525 

-.71 

-8.250 

-112.8 

520 

25 

2.4 

-.  1005 

-.22 

-16.219 

— 4l.O 

849 

26 

2.5 

-.  0962 

.29 

-16.060 

42.6 

755 

27 

2.6 

-.0503 

.58 

-8.637 

98.6 

324 

28 

2.7 

.0100 

.57 

1.945 

io4.o 

-211 

29 

2.8 

.0558 

.31 

10.555 

61.6 

-595 

30 

2.9 

.0693 

-.04 

13.468 

-4.3 

—666 

31 

3.0 

.0501 

-.  32 

9-997 

-61.1 

-426 

t©1 

=  329048.34  z©2=  6660.992256  £©©  =  -0.428921 

2©©  =  190.3603  E©@  =  -82  .  75361  2@©  =  16887.46 

E0©  =  4691.967  E©@  =  +98.2381628  Z  02=  2.66225444 

2©©  =  83.67703  z©©  =  -334446.428  £©©  =  -4835.7037 

£©©  =  23599.0  z@2=  +159.8717 

■- 

1 

1 
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TABLE  VI.-  TAYLOR*  S  SERIES  ITERATION  APPLIED  TO  EXAMPLE  III 


3.1852 

-8.5878 


-10.8324 


-10.8027 

-17.4231 

-15.2285 

-6.4301 

6.5720 


13.2823 


-13.8916 

-21.3879 


-10.1541  -18.2064 
-19.3557  -7.3413 


=  4286. 39139342 


2  (8X9)  =  -228.77281991 


£  (A;  (10)  =  35.627944979 
z  ®(n)  =  0.754833829 
2  ®  ©  =  1.43276343 
£  (9)  2  =  4l4o. 34148937 

£  (9) (10)  =  -  4.987891922 


£  @@  =-35.532782573 

z  @©  =  51.02915144 


£  Ud)a  =  0.3725473211 


Lg)  @  =  0.0542286301 


£  fid)  0.3)  =  -0.014427701 


£  @2  =  0.3777271561 


20®  =  -0.438887876 


tadjuk  V 11 . —  IRONY’S  METHOD  APPLIED  TO  EXAMPLE 
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O  .2  .4  .6  .8  1.0  1.2  1.4  1.6  1.8  2.0 

Independent  variable,  t 

(b)  Response  qm(t)  to  the  pulse  F(t) . 


Figure  /.-  Pulse  input  and  associated  response  used  in  example  I. 


Figure  3.  -  Response  to  a  step  input  used  in  example 
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Figure  5- First  derivative,  q  (t),  of  the  response,  q  (t),  used  in  example 
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Figure  6 -  The  variation  with  time  of  four  quantities  used  in  example  JH . 


Abstract 


A  least-squares  method  for  calculating  coeffi- 
°f  f  1inear  differential  equation  directly 
from  transient-response  data  is  presented.  Examples 
illustrating  the  application  of  the  method  to  the 

thP  «!  ?°n  °  aircraf tab ili ty  parameters  from 
the  airplane  response  to  an  elevator  deflection  are 
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